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Abstract. In this paper, we derive Darboux transformation of the inhomogeneous Hirota 
and the Maxwen-Bloch(IH-MB) equations which is governed by femtosecond pulse propagation 
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tion is used to derive sohton solutions, positon solutions of the IH-MB equations. 
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In recent years, nonlinear science has emerged as a powerful subject for explaining the mys- 
tery present in the challenges of science and technology today. Among nonlinear science, the 
interplay between dispersion and nonlinearity gives rise to several important phenomena in 
optical fibers, including parametric amplification, wavelength conversion, modulational insta- 
bility (MI), soliton propagation and so on. Among all concepts, solitons, positons and rogons 
have been not only the subject of intensive research in oceanography [T|[2] but also it has 
been studied extensively in several areas, such as Bose-Einstein condensate, plasma, superfluid, 
finance, optics and so on [3j-[9j. 

An important ingredient in the development of the theory of soliton and of complete inte- 
grability has been the interplay between mathematics and physics. In 1973, Hasegawa and 
Tappert [lOj modeled the propagation of coherent optical pulses in optical fibres by nonlinear 

^Corresponding author. 
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1. Introduction 
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Schrodinger (NLS) equation without the inclusion of fibre loss. They showed theoretically that 
generation and propagation of shape-preserving pulses called solitons in optical fibres is possible 
by balancing the dispersion and nonlinearity. 

In 1967, McCall and Hahn yjj had explained a special type of lossless pulse propagation 
in two-level resonant media. They have discovered the self-induced transparency(SIT) which 
can be explained by Maxwell- Bio ch( MB) equations. The coherent absorption takes place and 
the media becomes optically transparent to that particular wavelength when the energy dif- 
ference between the two levels of the media coincide with the optical wavelength. Burtsev 
and Gabitov |12| have considered MB equations with pumping and damping which is useful in 
optical pumping during the propagation of optical pulses in resonant atoms, and in their paper 
the Lax pair was presented for deformed MB systems. 

The important constraint to the NLS soliton namely the optical losses can be somewhat 
compensated with the effect of SIT. Therefore the system will be governed by the coupled system 
of the NLS equation and the MB equation (NLS-MB equations) if we consider these effects for a 
large width pulse. The coexistence of NLS solitons and SIT solitons in erbium-doped resonant 
fibres was experimentally observed by Nakazawa et al 13,14]. Recently, multi-soliton solutions 



of coupled NLS-MB equation was shown in |15|. In 16 , the periodic solutions have been 



generated through Darboux transformation and later rogue wave solutions were derived from 



breather solution in 17 . Modeling photonic crystal fiber for efficient soliton pulse propagation 



at 850 nm was surveyed in 18 . It presents new types of Dark-in-the-Bright solutions also called 
dipole soliton for the higher order nonlinear NLS (HNLS) equation with non-Kerr nonlinearity 
under some parametric conditions and subject to constraint relation among the parameters 
in optical context in [l9]. Impact of fourth-order dispersion in the modulational instability 
spectra of wave propagation in glass fibers with saturable nonlinearity was considered in [20| . 
The HNLS-MB equations as a higher-order correction of NLS-MB equations were shown that 
they allow soliton-type pulse propagation under a particular parametric condition [21,22 . 



For a reduced dynamical equation, the erbium-doped fibre system was proven to allow soliton- 
type pulse propagation with pumping [23j. The Lax pair and the exact soliton solution for 
Higher-order nonlinear Schrodinger and Maxwell-Bloch(HNLS-MB) equations with pumping 
was derived in |23| . 

Kodama f24] has shown that with suitable transformation and omitting the higher-order 
terms, higher order Nonlinear Schrodinger equation equation can be reduced to the Hirota 



equation 25 whose rogue wave solution is already reported in 26,27 . In a similar way,after 



suitable choice of self steepening and self frequency effects, the HNLS-MB equations can be 
reduced to a coupled system of the Hirota equation and MB equation [28j. The H-MB equations 
can be seen as the higher order correction of the NLS-MB equations and is the coupled system 
of the Hirota equation and the MB equation |31|. The H-MB system has been shown to be 
integrable and also admits the Lax pair and other required properties for complete integrability 

It is well known that the Darboux transformation is an efficient method to generate the 
soliton solutions for integrable equations [29||3^. The determinant representation of n-fold 
Darboux transformation of AKNS system was given in j^. In it constructed n- folds 
Darboux transformation of the H-MB equations, meanwhile the rogue wave solutions of the 
H-MB equations were obtained using the Darboux transformation. 

Inhomogeneous integrable equations become more and more attractive |^. Recently, K. 
Porsezian and C. G. Latchio Tiofack, Thierry B. Ekogo, etal. consider dynamics of bright soli- 
tons and their collisions for the inhomogeneous coupled nonlinear Schrodinger-Maxwell-Bloch 
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equations which describes propagation of an optical sohton in an inhomogeneous nonhnear 
waveguide doped with two level resonant atoms |33|. Soliton interactions in a generalized in- 
homogeneous Hirota-Maxwell-Bloch(IH-MB) system were considered in |34[|38| with symbolic 
computation but positon solutions of IH-MB is still unknown. 

The purpose of this paper is to derive the determinant representation of Darboux transfor- 
mation which is used to derive soliton solutions, positon solutions of the IH-MB equations. 

The paper is organized as follows. In Section 2, the Lax representation of IH-MB equations 
will be introduced firstly. In Section 3, we derived the one-fold Darboux transformation of 
the H-MB equations. In Section 4, the determinant-formed generalization of one-fold Darboux 
transformation to 2-fold Darboux transformation of the IH-MB equations will be given. Using 
these Darboux transformations, one soliton, two soliton and positons are derived in Section 
5 and Section 6 by assuming trivial seed solutions. Section 7 is devoted to conclusion and 
discussions. 



2. Lax representation of the IH-MB system 
In this paper, we will concentrate on the inhomogeneous Hirota and the Maxwell-Bloch(H- 



MB) system as following specific form 34,38 



= -{ai{z)Et + a2{z)E + ias{z)Eu + a^{z)Eut + a,{z)\E\^Et (2.1) 
+iaQ{z)\E\^E + a7{z)p), 
Pt = 2bi{z)Er]-2ib2{z)ujp, (2.2) 
Tit = -h^{z){Ep^ + E''p), (2.3) 

with constraint 

dhi 

^2 = as = 6a46?; = 2a^hl (2.4) 

In the equations above, z and t represent the normalized distance and time respectively, E{z^ t) 
denotes the slowly varying envelope axial field, p{z^ t) is the measure of the polarization of the 
resonant medium, and r]{z^t) represents the extent of the population inversion. ai{z) results 
from the group velocity and a2{z) describes the amplification or absorption. The coefficients 
as{z) ^ aQ{z) represent the group velocity dispersion (GVD), the third-order dispersion(TOD) 
[35] , self-steepening (SS) [36j, and self-phase modulation respectively, a^^z) is the parameter 
describing the averaging with respect to inhomogeneous broadening of the resonant frequency. 
bi{z) and 62 (^) depict the character of interactions between the propagation field and the erbium 
atoms. The real parameter is a constant corresponding to the frequency, and the * denotes 
the complex conjugate. If we set 

ai = a2 = 0, as = -^^^ ^4 = = -6/3, Qq = -a, = -2, bi = 1, 62 = -1, (2.5) 

the inhomogeneous H-MB equation will be reduced to H-MB equation as following 



Pt 
Vt 



1 



\E\'E) 

2iujp + 2Er}, 
-{Ep* + E*p). 



/3{Eut + Q\E\^Et) + 2p, 



(2.6) 

(2.7) 
(2.8) 



We will call the inhomogeneous Hirota and the Maxwell-Bloch system when a = 2,/3 = — 1 
the classical H-MB equation. The linear eigenvalue problem of IH-MB takes the form 



= U^, (2.9) 
= V^, (2.10) 

where U and V can be expressed in following polynomials about complex constant eigenvalue 
parameter A 

- a(J + '^^f)-X...U„, (2.11) 

_ X 3 /-4a4 0\ ,2/ -2^3 -4a4biE\ ,/ -ai - 2a4bl\E\'^ -2ihi{a^E - ia^EtfX 
4a4y \AaihiE* 2ia3 ) ^ \2ia:ihiE* - 2aihiE; ax^2aihl\E\^ ) 

= AV3 + AV2 + XV1 + V0+ ^ V V.i, (2.12) 
where 



2 \-p -''7 

^1 = -^vVT^^' ^1 •= '^o.AE*E^ + aiE + ia3^t + a^Eu, (2.14) 

Ci := iasblEE'' + a^bKE'Et - EE;). (2.15) 
\^ denotes the coefficient matrix of term and 

* = *(A)=f!:fM'!l) (2.16) 



$2(A,t,z) 

is an eigenfunction associated with eigenvalue parameter A of linear system eq.(2.9 ZTo] ). 

Using the linear equations of H-MB equations, One- fold Daroux transformation for IH-MB 
equation will be introduced in the next section. 

3. One-fold Daroux transformation for the IH-MB equation 

In this section, we will give the detailed proof of the one-fold Daroux transformation for the 
IH-MB equation. Firstly, we consider the transformation about linear function $ 

^/ _ r$ = (A^-5)$, (3.1) 

where 

A = ""A , S = 'A . (3.2) 

V«21 «22 / V^21 522/ 

New function is supposed to satisfy 

= (3.3) 

= V'^'. (3.4) 
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Then matrix T can be proven to satisfy following identities 

Tt + TU = U'T, 
n + TV = V'T. 



(3.5) 
(3.6) 



Bring the form of matrices A and S into eq.(3.5) and comparing the coefficients of both sides 
will lead to following condition 

ai2 = a2i = 0, (aii)t = (022)* = 0. (3.7) 

Therefore we will choose A = I and T = XI — S in the following part of this paper. The relation 
between E,p, rj and new solutions E',p', rj' which is called Darboux transformation can be got 



by eq. (3.5) and eq. (3.6). 



From (3.5), we have 
E' 



E + 2b^hi2, 

hE 
biE* 



S-S 







biE 



-biE* 



[S, aslS, 



and S should have a condition as S21 = s^2- By (3.6), following identity can be got 

1 



-S, + {XI - 5)(AV3 + AV2 + Al^i + Vo + 



X + iub2 



V-i) 



{x^V3 + AV2' + xv; + + 



1 



A + iub2 



VL,)iXI-S). 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



Multiplying both sides of eq.(3.10) hj XI — S can lead to 

-S,{X + iub2) + {XI - 5)(A3(A + lub2)V3 + X\X + iub2)V2 + A(A + iub2)Vi + Vo{X + iub2) + V.i] 
= (A3(A + iub2)V3 + X\X + iub2)Vi + A(A + iub2)Vi + V^{X + iub2) + VLi){XI - S). 
For term with A°, we get following identity 

S,iujb2 + Siiujb^Vo + V-i) = {iujb2V^ + VLi)S, 

which further leads to 

s, = {y^-iuj-\^vLi)S - s{-iuj-\^v.i + Vo). 

For term with A, we get following identity 

S, = {iub2Vo + V.i)- S{iub2Vi + Vo) + {iujb2Vi + V^)S+{-iub2V^-VLr). 
For term with A^, we get following identity 

{iub2Vi + Vo) - S{iub2V2 + Vi) = {iub2V[ + yj) - {iojb2Vi + V[)S. 
For term with A^, we get following identity 

{iujb2V2 + 1^1) - S{iujb2V^ + V2) = {iujb2V^ + Vl) - {iujb2V3 + V^)S. 

For term with A^, we get following identity 

V^ = V2-[S,V3], 

From above several identities, we can get 



E' = E + 2b^hi2, 
5 



(3.12) 



VL^ = (5 + iujh2)V-i{S + iujh2)~^, (3.13) 

which gives one fold transformation of one-fold Darboux transformation of H-MB equations. 
Suppose 

S = HAH-^ 



(3.14) 



where A 



Ai 






A2 



H 



^i{Xi, t,z) ^i{X2, t,z) 
$2(Ai,t,;^) $2(A2,t,^) 



$1,1 

$2,1 



$1,2 
$2,2 



12) 



H 



(3.15) 
(3.16) 



In order to satisfy the constraints of S and make V_^i having similar form as i.e. ^21 
foUowing constraint wiU be considered 

^$i(Ai,t,z) -$*(Ai,t,z) 
$2(Ai,t,z) ^l{Xut,z) 

The detailed determinant form of one-fold Darboux transformation of IH-MB equations in 
form of eigenfunctions will be given in the next section. 

4. Determinant representation of Darboux transformation 

In this section, we will give determinant representation of the first two Darboux transfor- 
mation of the IH-MB equations. Other higher-order Darboux transformation can be got in 
similar way which can be seen clearly in our paper ^31j. Firstly, we introduce n eigenfunctions 
$1 A 

'M = $(A = A^), i = 1, 2 with constraints on eigenvalues as A2i-i = —Ki the reduction 



conditions on eigenfunctions as $2,2^ = $1, 



2i-l5 



2,2i-l 



2r 



As the simplest Darboux transformation, the determinant representation of one-fold Darboux 
transformation of the IH-MB equations will be given in the following theorem. 

Theorem 4.1. The one-fold Darboux transformation of the IH-MB equations is as following 



where 



Ti(A,Ai,A2) = A/ + t. 



[1] 



1 

a; 



(Ti)n 

(Ti)21 



(Ti)i2 
(Ti)22 



(4.1) 



and 



( 



411 - 



1 

a; 



$2,1 
$2,2 

$2,1 
$2,2 



Al$l,l 
A2$l,2 

Al$2,l 
A2$2,2 



$1,1 
$1,2 

$1,1 
$1,2 



Al$l,l 
A2$l,2 

Al$2,l 
A2$2,2 



Ai = 



$1,1 
$1,2 



$2,1 
$2,2 



(4.2) 







1 





A 







1 







(Ti)ii 




$1,1 


$2,1 


Ai$i,i 


, (Tl)l2 = 


$1,1 


$2,1 


Ai$i,i 


(4.3) 






$1,2 


$2,2 


A2$l,2 




$1,2 


$2,2 


A2$l,2 








1 













1 


A 




(Tl)21 




$1,1 


$2,1 


Al$2,l 


, (Tl)22 = 


$1,1 


$2,1 


Al$2,l 


(4.4) 






$1,2 


$2,2 


A2$2,2 




$1,2 


$2,2 


A2$2,2 




f/[i] 


= V 


+ [Ti 






^1 A=-iw62^- 




\=—i(joh2 ) 


(4.5) 



= E + 2b^hi2 ^ E - 2b^ 



1 m) 
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Ai ' 



[1] _ 



2?7(ri)ii(Ti)i2 -p*(Ti)i2(ri)i2 +p(ri)n(Ti )n 
(Ti)n(ri)22 - (Ti)i2(ri)2i 



\X=—ioub2 ) 



(4.6) 
(4.7) 



[1] ^ n{{Ti)ii{T,)22 + (Ti)i2(ri)2i) - p*(ri)i2(Ti)22 +p(ri)ii(Ti)2i 
^ (Ti)n(ri)22 - (Ti)i2(ri)2i ''=-''^''^- ^ ^ 

We can find the transformation Ti iias following property 



Ti(A;Ai,A2)| 



\=\i 



^1, 

$2, 



= 0, 



(4.9) 



where i = 1, 2. 

This one-fold transformation will be used to generate one-soliton solution from trivial seed 
solution of the IH-MB equation. 

In the next part, we will generalize the Darboux transformation to two- fold case which is 
contained in the following theorem. 



Theorem 4.2. The two-fold Darboux transformation of H-MB equation is as following 



where 



T2{X; Ai, A2, A3, A4) = X'l + t?X + C = 



J21 



1 AT2)i1 (T2)l2 



(T2) 



21 



T 



2)22 



(4.10) 



(T2)ll = 



Ao = 



$1,1 $2,1 Ai$i,i Ai$2,i 

$1,2 $2,2 A2$l,2 A2$2,2 

$1,3 $2,3 A3$i,3 A3$2,3 

$1,4 $2,4 A4$i,4 A4$2,4 



(T2)21 = 



We can find 



^2(A; Ai, A2, A3, A4)| 



$1,^ 

$2,i 







(4.11) 





1 





A 





A2 







1 





A 







$1,1 


$2,1 


Al$l,l 


Al$2,l 


A?$i,i 




$1,1 


$2,1 


Al$l,l 


Al$2,l 


A?$l,l 




$1,2 


$2,2 


A2$l,2 


A2$2,2 


Ai$i,2 


, (T2)l2 = 


$1,2 


$2,2 


A2$l,2 


A2$2,2 


Ai$i,2 




$1,3 


$2,3 


A3$l,3 


A3$2,3 


Ai$i,3 




$1,3 


$2,3 


A3$l,3 


A3$2,3 


Ai$i,3 




$1,4 


$2,4 


A4$l,4 


A4$2,4 


Al$l,4 




$1,4 


$2,4 


A4$l,4 


A4$2,4 


Al$l,4 




1 





A 













1 





A 


A2 




$1,1 


$2,1 


Al$l,l 


Al$2,l 


A?$2,l 




$1,1 


$2,1 


Al$l,l 


Al$2,l 


A?$2,l 




$1,2 


$2,2 


A2$l,2 


A2$2,2 


Ai$2,2 


, (T2)22 = 


$1,2 


$2,2 


A2$l,2 


A2$2,2 


Ai$2,2 




$1,3 


$2,3 


A3$l,3 


A3$2,3 


Ai$2,3 




$1,3 


$2,3 


A3$l,3 


A3$2,3 


Ai$2,3 




$1,4 


$2,4 


A4$l,4 


A4$2,4 


A4$2,4 




$1,4 


$2,4 


A4$l,4 


A4$2,4 


Al$2,4 



(4.12) 



where i = 1,2,3,4. Similarly, for transformation 

T2t + T2U -- 
T2, + T2V -- 



T2, following transformation formula holds 

= [/[^ITs, (4.13) 
= yP]T2, (4.14) 



by which the relation between E^p^rj and E^'^\ p^'^\ rj^'^^ wiU be got in the foUowing relation 



[2] 



(4.15) 



T2\\=-iujh2^-iT2 ^\\=-iujh2- (4.16) 

This gives the relation between E^p^r] and E"'^!,])'^!, ry'^l One can also get following two-fold 
Darboux transformation in detail. 



[2] 



V 



[2] 



^-((tr')i2, 

277(T2)ii(r2)i2 -p*(r2)i2(T2)i2 +p(T2)ii(r2 )ii 

(^2)11(^2)22 — (^2)12(^2)21 

'7((r2)ii(r2)22 + (r2)i2(r2)2i) -p*(r2)i2(r2)22 +P(r2)n(r2)2i 



\X=—iu:b2 5 



(4.17) 
(4.18) 

(^2)11(^2)22 - (r2)i2(T2)2i ''=-^'^''^' ^^-^^^ 

where (tf )i2 is the element at the first row and second column in the matrix of . This 
transformation will be used to generate two-soliton solutions of the IH-MB equation later. 

As an application of these determinant representation of Darboux transformations of IH-MB 
equations, soliton solutions and positon solutions will be constructed in the next section. 

5. Soliton solutions of the IH-MB equations 

In this section, first, we will consider the construction of one soliton solution of the IH-MB 
equations with suitable seed solutions. Bring trivial seed solutions as = 0,p = 0, 77 = 1 into 



linear equations eqs.(2.9 2.10)., then the linear equations become 



where 



$ = 
U = 
V = 



$1 
$2 

A 

-A, 

— 4A^a4 — 2iA^a3 — Aai 




bia7 







hia7 



) 



4A^a4 + 2iX'a3 + Aai + 2X^3^^/ 
Easy calculation can lead to following eigenfunctions 

$1 = exp(At + / -4A^a4 - 2iA^a3 - Aai - 



$9 



exp(-At + / 4A^a4 + 2iA^a3 + Aai + 



2A + 2iub2 

bitty 



ttz H ), 



dz — 



2 



+ 16), 



(5.1) 
(5.2) 

(5.3) 
(5.4) 
(5.5) 

(5.6) 
(5.7) 



2A + 2iujb2 " 2 

where XQ^y^ and 6 are all arbitrary fixed real constants. Substituting these two eigenfunctions 



into the one-fold Darboux transformation eq.(4.6), eq.(4.7) and eq.(4.8) and choosing A = 



tti + /5ii, a;o = 0, j/o = 0, ^ = 0, then the following solition solution are obtained: 



El = 



Pi 



2ai 



Az 



2{ai + j3ii + ujh2i) 



2toi 



-Az 



2{-ai + l^ii + ujh2i) 



) 



(/5i+cj62)^ 



at 



/ 2tol ^rT^^ , . 



-2tai- 



-2tai- 



y.l-\-Pli-\-Lob2i 



sech^{ 



Az 



2{ai + + cl;620 



2tai 



-Az 



2{-ai + + u)b2i) 



)), 



r]i = {{cosh?{- 



Az 



sech 



2(«i + Pii + ujb2i) 
Az 



2tai 



-Az 



^2(q;i + + ujh2i) 



- 2toi + 



2(-q;i + l^ii + c<j52^) 
2(-«i + /3ii + ub2i) 



) 



at 



(/5i+a;62)2 + «?) 



) 



where 



= 8a4ai + 2miaiCc;62 — 32^4^1/3^ — 48a4a^/3^ — 24:a4^a\PiUjh2 — 12^3^1/5^ 

+32m4a;L/5i + 8a4/3;L + 8a4/3;L^^2 + 4:iasai — 12asail3i — 4:asaiUjb2 — 24^40^1/3^(^62 
+8ia4al(jjb2 + 4a3/3^ + 4:asf3iUjb2 + 2aia\ + Aiaiai^i — 8iasai(3i(jjb2 
—2ai0l ~ 2aiPiUjb2 + ftiay). 
Similarly, substituting these two eigenfunctions into the one-fold Darboux transformation 



eq.(4.6), eq.(4.7) and eq.(4.8), and taking ai = z,a3 = — l,a4 = z,a7 = = 1,62 = = 



1.5,^1 = 0.5, /3i = 2, then the one-solition solutions of the classical H-MB equations can be 
obtained whose evolution is given in Figjl} which clearly indicates that E and p are bright 
solitons because their waves are under the flat non- vanishing plane whereas 77 is a dark soliton. 




Figure 1. One solition solution (E^p^rj) of the IH-MB equations when ai = 
z, as —1, a4 = z, a7 = z, 61 = 1, 62 = ^, = 1.5, ai 0.5, /3i 2,. 



Now let us discuss about the construction of the two-soliton solution of IH-MB system. For 

the purpose of construction of two soliton solution, we need to use two spectral parameters 

Xi = ai + /3ii and A2 = 0^2 + /32i- After the second Darboux transformation, we can construct 

the two solition solution. As the general form of two soliton solution is tedious in nature, for 

simplicity, we will give only the two soliton solution of E taking values as ai = z^a^ = — 1, a4 = 

z, a7 = z, 61 = 1, 62 = ^, ^ = 1-5, ai 0.5, /3i = 2, a2 = 1, 1^2 1-5 in appendix. 
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We also construct two soliton solutions for p and ?] in a similar manner. For completeness, 
instead of giving complicated forms of p and rj, the graphical representation of them is shown 
in Figj2} 




(i^n f mmmi 



{\p? 




iv) 




Figure 2. Two solition solution (E^p^rj) of the IH-MB equations when ai = 
aa — l,a4 = z^a^ = z^bi = l^b2 = z^uj = 1.5,^1 0.5, = 2,^2 = '^^^2 
= 1.5. 

If we suppose ai = 0, a2 = 0, = —1, = 1, as = 6, = —2, = —2, bi = 1, 62 = —1, this 



case wiU go to the classical H-MB equation 31 with constant coefficients. 



6. Bright and dark positon solutions of IH-MB system 

For the two soliton solution constructed in the last section, if the second spectral parameter 
A2 is close to the first spectral parameter Ai, doing the Taylor expansion of wave function up 
to first order near Ai will lead to a new kind of solution. This is exactly the so-called positon 
solutions. In this section, the construction of positon solution of IH-MB equations will be given. 
Firstly, following four linear functions out of linear system will be used to construct the second 
Darboux transformation which further generate the positon solutions. 



$1 
$2 

$3 
$4 



= exp(Ait + 
exp(— Alt + 
exp(A2t + 

= exp(-A2t + 



4A^a4 — 2iXlas — Xiai — 
4:Xla4 + 2iXlas + Aiai + 
4A2a4 — 2iA2a3 — A2ai — 
4A2a4 + 2zA2a3 + A2ai + 



bia^ 



^2,3 



2Ai 


+ 2iub2 




bitty 


2Ai 


+ 2iub2 




bitty 


2A2 


+ 2iub2 




bitty 


2A2 


+ 2iujb2 


■=4 


>2; $2,2 : 


■=4 


>4; $2,4 : 



dz + 
dz — 
dz + 

dz — 



2 ' 
xq + iyo 
2 

xq + iyo 
2 

xo + iyo 



(6.1) 
(6.2) 
(6.3) 
(6.4) 

(6.5) 
(6.6) 



Then we define the following functions as 

$1,1 := 01; $1,2 := 
$1,3 := 03; $1,4 := 

Now we take A2 = Ai + e{l + i) and using the Taylor expansion of wave function ^3 and 04 up 
to first order of e in terms of Ai. Substitution of these manipulations into the second Darboux 
transformation discussed in the last section will lead to positon solutions. For example, after 
taking values as ai 0.5, a3 = z,a4 = z^a^ = 1, 61 = 1, 62 = 1, ^ = 1.5,^1 0.5, ;5i 1, 
i.e. the case of classical H-MB equations |31|, the positon solutions (Ep^pp^rjp) can be derived. 
Here for simplicity, we only give positon solution Ep in following form 
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Ep = -0.2e2^^+2-^^^'"°-^^^^^^^^^^^^[(-130000000000z^ + 140000000000iz^ + 11420118340z 
+591715977iz - 20000000000t - 20000000000)e-'-^-™™°^'+°-^^^^2^°^^^^ 
+ (591715977iz + 140000000000iz^ + 130000000000^^ - 11420118340z + 20000000000t 
— 20000000000)e^^'''^^^^^^^^^^^~^'^''^^^^^^^^^] 

/(2000000000cos/i(-15z2 + 1.153846154z - 2t) + 4000000000^^ + 52000000010tz^ 
+2000000000 - 4568047338tz + 365000000000^^ 
-28035502960^^ + 1307692307^^). 
The pictorial representation of positon solutions of the IH-MB equations is shown in Figj3} 




Figure 3. One positon solution [E^p^r]) of the IH-MB equations when ai = 
0.5, a3 = a4 = a7 = 1, 6i = 1, 62 = 1, = 1.5, ai 0.5, /3i 1. 

whose density plot is as Figji} Next, after taking values as ai = 0.5, a3 = e^,a4 = e^,a7 = 




Figure 4. One positon solution (E^p^rj) of the IH-MB equations when ai = 
0.5, as = z^a^ = z, = 1, 61 = 1, 62 = 1, ^ = 1-5, 0^1 0.5, /3i 1. 



1,61 = 1,62 = 1,^ = 1.5,^1 = 0.5, ;5i = l,the picture of positon solutions of the IH-MB 
equations is as FigjSj whose density plot is as FigjH) 

From above we find that E and p are bright positon solutions whereas 77 is a dark positon. 
In a similar way, using the higher order Darboux transformation, one can also generate higher- 
order bright and dark positon solutions which will be omitted here. These positons with variable 
coefficients are different from the classical H-MB equations which can be seen from their graphs. 
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Figure 5. One positon solution [E^p^r]) of the IH-MB equations when ai = 
0.5, a3 e^, a4 e^, a7 = 1, 61 = 1, 62 = 1, ^ = 1-5, oii 0.5, /3i 1. 





Figure 6. One positon solution (E^p^rj) of the IH-MB equations when ai 
0.5, as e^, e^, = 1, 61 = 1, 62 = 1, = 1.5, o^i 0.5, /3i 1. 



7. Conclusion and Discussions 



In this paper, we derived the Darboux transformation of the inhomogeneous Hirota and the 
Maxwell-Bloch(IH-MB) equations governed by ultra-short pulse propagation through erbium 
doped optical waveguide. Further matrix representation of Darboux transformation of this 
system is constructed. As examples, soliton solutions, positon solutions of the IH-MB equations 
have been constructed explicitly by using Darboux transformation from trivial solutions seed 
solutions. There are a few unclear interesting questions such as the physical interpretations 
and observation of higher-order positon solutions, rogue waves solutions and their applications 
in physics? 
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8. Appendix 



^2-sol 

= [8(17 + 24z + 9z')Hl3 + 18z + Q^^)l-¥e''-'->'^^^^^ 

_^8(13 + 18^ + 9^2)§-fv-8^+^+3^^+§^^^^^^(4+3z)+23z^-f iz-6iz^/ l + 4:i + 3iz 2 

Vl7 + 24z + 9z2^ 

_2(13 + 18z + 9^2)1(17 + 24z + Q^^)h-I^e^^^-''^^^^^^^^^ 

Vl3 + 18z + 9z2' 

-2(17 + 24z + 9x;^)§-|v-i^^+^^+4^^+f^^^^^^(i+f^)+^^^^-f 2 + 3i + 3iz 4 

V13 + 18z + 9z2^ 

+4i(17 + 24z + 9z')Hl3 + 18z + 9^^)l-¥e'^-'^'''->^^^^ 

^ ^ ^ ^ Vl7 + 24z + 9z2^ 

+4zfl7 + 24Z + 9Z^)^- - l,->-12^+2mit+|arctan(| + |z)+22z^-f i^+2fe^/ 2 + 3z + 3lZ 4 

^ ^99^ Vl3 + 18z + 9z2^ 

— -\- ISz + 9Z^)^- — l^^g-8^+*+3^^+|c^^ctan(4+32;)+232;2-13/3iz-6i2;2 / 1 + 4i + 3iz 



2 



9 3 V17 + 24^ + 9^2^ 

-4i(13 + 18^ + 9^2)i(17 + 24Z + 9^2)|-|igl2.-2mit-farcton(| + |.)-22.^-f ^.+24».^( 2 + 3i + 3iz 

V 13 + 18;^ + 92;2 

2 + 3i + 3iz ,4, 1 + 4:i + 3iz x2, „ ^ r>^2,^„ „ 9n i 

/ ( , = 9 -8 13 + 18^ + 9^2)9 17 + 24:z + 9z^)^cos 

V13 + 18^ + 9^2^ V17 + 24^ + 9^2^ 

4 3 3 1 2 4 

{t + -arctan{- + -z) + -ln{13 + 18z + 9z'^) - lOz + 30z'^ - -arctan{A + 3^) - -/n(17 + 24^ + 9^^)) 
9 2 2 3 9 9 

_^5(^]_7 _^ 24z + 9^2^|g-42;+t+|arcfan(| + |2;)-|arctan(4+32;)-2;2 _^ ^_202;+3t+|arctan(| + |2;) + |arctan(4+32;)+452;2 
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_^5(13 + 18z + 9^2^j|g42:-t-|arctan(| + |2;)+|arctan(4+32;)+2;2 

+(13 + 18z + 9z')t (17 + 24z + Qz^)le^^'-^'-h--'--(l+l^^^^^ 
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